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: ABSTRACT 

I Sulfur gases are common components in the volcanic and biological emission 

I I on Earth, and are expected to be important input gases for atmospheres on 

^ I terrestrial exoplanets. We study the atmospheric composition and the spectra of 

^ I terrestrial exoplanets with sulfur compounds (i.e., H2S and SO2) emitted from 

their surfaces. We use a comprehensive one-dimensional photochemistry model 
and radiative transfer model to investigate the sulfur chemistry in atmospheres 
ranging from reducing to oxidizing. The most important finding is that both 



O I H2S and SO2 are chemically short-lived in virtually all types of atmospheres on 

^ I terrestrial exoplanets, based on models of H2, N2, and CO2 atmospheres. This 

^ I implies that direct detection of surface sulfur emission is unlikely, as their surface 

^ I emission rates need to be extremely high (> 1000 times Earth's volcanic sulfur 

emission) for these gases to build up to a detectable level. We also find that 
sulfur compounds emitted from the surface lead to photochemical formation of 
^ I elemental sulfur and sulfuric acid in the atmosphere, which would condense to 

form aerosols if saturated. For terrestrial exoplanets in the habitable zone of Sun- 
like stars or M stars. Earth-like sulfur emission rates result in optically thick haze 
composed of elemental sulfur in reducing H2-dominated atmospheres for a wide 
range of particle diameters (0.1 - 1 /im), which is assumed as a free parameter 
in our simulations. In oxidized atmospheres composed of N2 and CO2, optically 
thick haze, composed of elemental sulfur aerosols (Sg) or sulfuric acid aerosols 
(H2SO4), will form if the surface sulfur emission is 2 orders of magnitude more 
than the volcanic sulfur emission of Earth. Although direct detection of H2S and 
SO2 by their spectral features is unlikely, their emission might be inferred by 
observing aerosol-related features in reflected light with future generation space 
telescopes. 
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Subject headings: radiative transfer — atmospheric effects — planetary systems 
— techniques: spectroscopic — astrobiology 

1. Introduction 

A large number of super Earths have been detected by radial velocity surveys and 
transit surveys. Super Earths are exoplanets with masses no more than 10 times the mass 
of Earth. The atmospheres of super Earths are important because characterization of super 
Earth atmospheres is a substantial step towards eventually characterizing truly Earth-like 
exoplanets. Attempts to observe super Earth atmospheres are growing (e.g., Batahla et al. 
2011 for Kepler 10 b; Demory ct al. 2012 and Ehrenreich et al. 2012 for 55 Cnc e), and 
one super Earth/mini Neptune GJ 1214b is being observed in as much detail as possible 
(e.g., Bean et al. 2010; CroU et al. 2011; Desert et al. 2011; Berta et al. 2012; De Mooij 
et al. 2012). In particular, the transmission spectra of GJ 1214b is nearly fiat from 0.6 
to 5 /im, which has ruled out a planet with an H2-dominated cloud-free atmosphere. The 
observational push to super Earth characterization has the potential to provide a handful 
of super Earth atmospheres to study in the coming years. In the more distant future, the 
community still holds hope that a direct-imaging space-based mission under the Terrestrial 
Planet Finder concept will come to existence, with planets and their atmospheres observed 
in reflected light. 

The concentration of trace gases in super Earth atmospheres is controlled by the com- 
ponent gas emission from the surface and subsequent sinks in the atmosphere (chemical 
reactions initiated by UV photolysis). Even with a trace amount, some gases may leave 
significant footprints in the planet's spectra, for instance H2O, GO2, and O3 in Earth's at- 
mosphere. In contrast to giant exoplanets where the atmospheric composition is mainly 
controlled by the elemental abundance, the steady-state composition of super Earth atmo- 
spheres is mainly controlled by photochemical processes. We have developed a comprehensive 
photochemistry model that computes chemical compositions of any terrestrial exoplanet at- 
mospheres, ranging from reducing to oxidizing, with all key non-equilibrium processes taken 
into consideration, including photolysis, chemical kinetics, vertical diffusion of molecules, 
atmospheric escape, dry and wet deposition, and condensation and sedimentation of con- 
cern condensable species (Hu et al. 2012; referred to as Paper 1 hereafter). The effects 
of surface gas emission on different super Earth atmospheres can be investigated with the 
photochemistry model. 

Sulfur gases emitted from the surface and their photochemical products significantly 
shape the spectra of rocky bodies in the Solar System. The most striking feature of Venus' 
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atmosphere is a high planetary albedo due to thick H2SO4 clouds. Photochemistry models 
of the Venusian atmosphere have been developed and the formation of H2SO4 in the dry 
C02-dominated atmosphere have been simulated (e.g., Yung & DeMore 1982; Zhang et 
al. 2012; Krasnopolsky 2012). The latest photochemistry model of Venus' atmosphere has 
assumed a constant mixing ratic0 of SO2 (i.e., ~ 10 ppm) at the altitude of 47 km (implying 
significant source of SO2 from below) and predicted the formation of H2SO4 at the altitudes 
around 66 km (Zhang et al. 2012; Krasnopolsky 2012). lo, the innermost moon of Jupiter, is 
believed to have very intensive and variable volcanic activity with SO2 emission, and the lo's 
atmosphere is dominated by photolysis of SO2 and subsequent formation and polymerization 
of elemental sulfur (e.g., Moses et al. 2002). On early Earth, the sulfur chemistry may have 
been very different from now, featuring the formation of sulfur aerosols as well as sulfate 
aerosols as the atmosphere was anoxic (Kasting et al. 1989; Pavlov & Kasting 2002; Kasting 
& Catling 2003; Zahnle et al. 2006). The formation of insoluble sulfur aerosols is believed to 
be critical for the record of mass independent fractionation that timed the rise of oxygen in 
Earth's atmosphere (e.g., Farquhar et al, 2000; Pavlov & Kasting 2002; Zahnle et al. 2006). 
In addition, organosulfur compounds, such as dimethyl sulfide (DMS) and methanethiol 
(CH3SH), have also been suggested to be biosignatures of the early Earth (Pilcher 2003). 
The greenhouse effect of SO2 has been suggested to have contributed to the warming of early 
Mars (e.g., Halevy et al. 2007), a proposition that has been challenged by photochemistry 
studies that predict sulfur and sulfate aerosol formation as a result of SO2 emission on early 
Mars and the anti-greenhouse effect of these aerosols (Tian et al. 2010). 

Terrestrial exoplanets could have a wide range of sulfur gas emission. Sulfur gases are 
common volcanic gases on Earth, and in some scenarios may be more prevalent on exoplanets. 
Sulfur is a major building block for rocky planets and the abundance of sulfur is about one 
sixteenth that of carbon in the Solar System (Lodders 2003). On present-day Earth, sulfur 
compounds, mainly in the form of HgS and SO2, are dominant volcanic gases in additional 
to hydrogen, oxygen and carbon emissions. In the Earth's magma, the dissolved sulfur has a 
weight percentage ranging from 10~^ to 10~^ (e.g., Wallace & Edmonds 2011) and degassing 
of sulfur compounds as the magma decompresses when rising to the surface provides a global 
volcanic sulfur flux of about 3 x 10^ S cm~^ s~^ (Seinfeld & Pandis 2006). 90% of current 
Earth's sulfur emission is in the form of SO2, whereas the H2S/SO2 ratios for individual 
volcanoes vary widely between 0.01 and 1 (see Holland 2002 and references therein). An 
intriguing fact is that the amounts of sulfur compounds in Earth's atmosphere is extremely 
low despite the substantial emission rates, due to very short chemical timescales. On Earth, 



^Mixing ratio is defined as the ratio of tire amount of a gas in a given volume to the total amount of all 
gaseous constituents in that volume. 
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the lifetime of H2S and SO2 in the troposphere is only 2 days, which makes the mixing ratio 
of these gases in the atmosphere very small (Seinfeld & Pandis 2006). In Earth's troposphere, 
the main sink of is the reaction with the hydroxyl radical OH (Lelieveld et al. 1997), 
and the main sink of SO2 is the removal by dry and wet deposition. 

Another reason that H2S and SO2 photochemistry is interesting is that Earth-based 
biological processes involve sulfur compounds. There arc multiple ways that life can produce 
H2S, including the reduction of sulfate (e.g., Watts 2000) and the disproportionation of sulfur 
compounds of intermediate oxidation states (e.g., Finster 2008). In general, 10 to 50 percent 
of the global H,S emission on Earth is as a metabolic byproduct, whereas over 70% percent 
of natural SO2 emission is volcanic, although on modern Earth 90% of atmospheric SO2 is 
anthropogenic (Seinfeld & Pandis 2006). Understanding the atmospheric response to H2S 
and SO2 emission will allow us to examine whether or not H2S is a potential biosignature 
gas on a terrestrial exoplanet. 

Previous investigation of sulfur photochemistry in the context of super Earth character- 
ization has been very limited. Sulfur compounds are generally not considered in most models 
of terrestrial exoplanet atmospheres. Des Marais et al. (2002) describe the spectral features 
of certain molecules such as H2O, N2O, O3 and CH4 in terrestrial exoplanets. Miller- Ricci et 

al. (2009) present spectra of super Earths under equilibrium chemistry with photochemistry 
estimates and without considering sulfur compounds. Zahnlc et al. (2009) investigate sulfur 
photochemistry in hot Jupitcrs and suggest that HS and S2 can be generated photochcm- 
ically from H2S. Kaltcncggcr & Sassclov (2010) study the effect of H2S and SO2 on the 
terrestrial planetary spectra and suggest that 1-10 ppm SO2 has potentially detectable spec- 
tral features that could indicate active volcanism. Kaltenegger & Sasselov (2010), however, 
do not link the surface emission of sulfur gases to the mixing ratio of sulfur species via photo- 
chemistry involving sulfur compounds, or consider formation of aerosols. Domagal- Goldman 
et al. (2011) study the chemistry of organic sulfur compounds that are strongly linked to 
biology (e.g., CH3SH) in anoxic N2 atmospheres and suggest that the most detectable feature 
of organic sulfur gases is their indirect photochemical product, ethane. Moreover, the fate 
of surface emission of sulfur compounds, e.g., H2S and SO2, is yet to be explored for H2- 
dominated atmospheres. H2-dominated atmospheres, like N2 or CO2 atmospheres with water 
vapor, could also maintain a habitable temperature at the surface through collision-induced 
absorption (Pierrehumbert & Gaidos 2011; Wordsworth 2012). 

In all, it is still largely unknown whether or not H2S and SO2 spectral features can be 
observed in the future in the atmosphere of an exoplanet and whether or not the emission 
rate of sulfur compounds on a terrestrial exoplanet can be inferred. In this paper we investi- 
gate the atmospheric chemistry resulting from H2S and SO2 surface emission in atmospheres 
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having very different oxidation states ranging from reducing to oxidizing. Here, "emission" 
means tlie mass ffux from tlie planetary surface to tlie atmospliere tliat may include volcanic 
sources and biogenic sources. We focus on terrestrial exoplanet atmospheres that include 
super Earths, since those planets hold the most interest on the path to finding and charac- 
terizing planets that might harbor life. In § 2 we briefly describe our photochemistry model 
and radiative transfer model, with a focus of treatments of sulfur chemistry and aerosols. 
In § 3 we describe the key chemistry pathways involving HgS and SOg in both reducing 
and oxidizing atmospheres. In § 4 we present the main results on sulfur photochemistry in 
terrestrial exoplanet atmospheres and the spectral features of HgS, SO2, and photochemical 
aerosols. In § 5 we discuss on whether or not HgS can be a biosignature gas on a planet with 
atmospheric conditions different from Earth's. We present our conclusions in § 6. 



2. Model 

2.1. Photochemistry Model 

We have developed a comprehensive photochemistry model to investigate atmospheres of 
terrestrial exoplanets, and validated the model by simulating the atmospheric compositions 
of current Earth and Mars (see Paper I). We now describe briefly the main features of the 
photochemistry model and the specifics for this work. 

The purpose of the photochemistry model is to compute the steady-state chemical com- 
position of an exoplanetary atmosphere. The system is described by a set of time-dependent 
continuity equations, one equation for species at each altitude. Each equation describes: 
chemical production; chemical loss; eddy diffusion and molecular diffusion (contributing 
to production or loss); sedimentation (for aerosols only); emission and dry deposition at 
the lower boundary; and diffusion-limited atmospheric escape for light species at the upper 
boundary. Starting from an arbitrary initial state, the system is numerically evolved to the 
steady state in which the number densities no longer change. Because the removal timescales 
of different species are very different, the implicit inverse Euler method is employed for the 
numerical time stepping. The generic model computes chemical and photochemical reactions 
among 111 molecules and aerosols made of O, H, N, C, S elements, and formation of sulfur 
(Sg) and sulfate (H2SO4) aerosols. The numerical code is designed to have the capacity of 
treating both reducing and oxidizing atmospheres. For the chemical and photochemical reac- 



tions, we use the reaction rates data from both the NIST database (http:/ /kinetics. nist.gov) 
and the JPL publication (Sander et al. 2011). We have also adopted relevant reaction rates 
from Kasting (1990), Yung & DeMore (1999), and Moses et al. (2002). Sulfur polymerization 
reaction rates still lack consistent experimental measurements, and we adopt the reaction 
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rates proposed by Kasting (1990) and those proposed by Moses et al. (2002). Both sets 
of sulfur polymerization reaction rates are speculative and they are widely discrepant, the 
effect of which will be discussed later in section 3.1. Ultraviolet and visible radiation in 
the atmosphere is computed by the (5-Eddington 2-stream method, with molecular absorp- 
tion, Rayleigh scattering and aerosol Mie scattering contributing to the opacity. For the 
stellar input spectrum we used the Air Mass Zero (AMO) reference spectrum produced by 
the American Society for Testing and Material^ for Sun-like stars, and used the simulated 
non-active M star spectrum from AUard et al. (1997) for quiet M stars. 

In this paper we use the photochemistry model to study the sulfur chemistry in atmo- 
spheres ranging from reducing to oxidizing on terrestrial exoplanets. We use H2 dominated 
atmospheres as the representative cases for reducing atmospheres, and we use Nj, and CO2 
dominated atmospheres as the representative cases for oxidized atmospheres that could be 
both reducing and oxidizing. With the photochemistry model, we simulated the chemical 
composition of Hg, Ng, and COg-dominated atmospheres, with sulfur compounds emitted 
from the surface at various rates. Parameters of the atmospheric models are tabulated in 
Table [T] and details of water, carbon, and oxygen chemistry have been described in Paper 
I. We will focus on the sulfur chemistry and photochemistry in terrestrial exoplanet atmo- 
spheres in this paper. 

We consider in our atmospheric chemistry models the O, H, and S bearing species and 
a subset of C bearing species. The gaseous molecules considered in this paper are H, H2, O, 
0(^D), O2, O3, OH, HO2, H2O, H2O2, CO2, CO, CH2O, CHO, C, CH, CH^, 'CR,, CH3, 
CH4, CH3O, CH4O, CHO2, CH2O2, CH3O2, CH4O2, C2, C2H, C2H2, C2H3, C2H4, C2H5, 
C2H6, C2HO, C2H2O, C2H3O, C2H4O, C2H5O, S, S2, S3, S4, SO, SO2, ^S02, ^S02, SO3, 
H2S, HS, HSO, HSO2, HSO3, H2SO4, and Sg, and the aerosols considered are Sg aerosols 
and H2SO4 aerosols. This set of species is comprised of common H, O, and C bearing species 
and photochemical products of H2S and SO2 emission. We assume a constant H2O relative 
humidity at the surface of 60% to mimic the supply of water vapor from a liquid water ocean. 
To reduce the stiffness of the system and improve the numerical stability, "fast" species with 
relatively short chemical loss timescales are computed directly from the photochemistry 
equilibrium. We consider in this work 0(^D), ^CH2, C2H, ^S02, and ^S02 as fast varying 
species. As such, the photochemistry model rigorously finds steady-state composition of 
the atmosphere starting with initial compositions without any sulfur compounds. Once the 
model converged to the steady state, we checked explicitly the mass conservation of O, H, 
C, N, S atoms and verified the choice of fast species to have been appropriate. We have 



^ http: //rredc. nrel.gov/solar/spectra/amO/ 
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Tablc 1: Basic parameters for the atmosphere models in this paper. The scenarios are H2, 
N2, CO2 dominated atmospheres on an Earth-sized atmosphere in the habitable zone of 
a Sun-like star. The planet has Earth-like volcanic emissions for hydrogen, oxygen, and 
carbon species, and various H2S and SO2 emissions. Note that we do not consider any 
biotic contribution to the dry deposition velocities, and the volcanic carbon emission is not 
proportionally increased with sulfur emission. 



Parameters 


Reducing Oxidized 


Oxidized 


Main component 


90%H2, 10%N2 


90%CO2, 10%N2 


Mean molecular mass 


4.6 28 


42.4 


Planetary parameters 






Stellar type 


G2V G2V 


G2V 


Semi-major axis 


1.6 AU 1.0 AU 


1.3 AU 


Mass 






Radis 






Temperature profile 






Surface temperature 


288 K 288 K 


288 K 


Surface pressure 


10^ Pa 10^ Pa 


10^ Pa 


Tropopause altitude 


120 km 13.4 km 


8.7 km 


Temperature above tropopause 


160 K 200 K 


175 K 


Maximum altitude 


440 km 86 km 


51 km 


Eddy diffusion coefficient 








6.3 X 10^ cm^ s~"^ 1.0 X 10^ cm^ 


6.8 X 10"* cm^ 


Minimum 


2.5 X 10* cm^ 3.9 X 10^ cm^ 8"^ 


2.7 X 10^ cm^ 


Altitude for the minimum 


107 km 17.0 km 


11.6 km 


Near the top of atmosphere 


7.1 X 10^ cm^ 1.1 X 10^ cm^ 


7.6 X 10* cm^ 


Water and rainout 






Liquid water ocean 


Yes Yes 


Yes 


Water vapor boundary condition 


/(HjO) = 0.01 /{HjO) = 0.01 


/(H^O) = 0.01 


Rainout rate* 


Eartli-likc Earth-like 


Earth-like 


Gas emission^ 








3 X 10^^ em^^ s^l 3 X lO^^ cm^^ s^^ 


N/A 


H2 


N/A 3 X 10^" cm"^ s^l 


3 X 10^" cm-^ s~^ 


CH^ 


3 X 10** em-2 s"^ 3 X 10* cm"^ 


3 X 10* cm-^ 




Vary Vary 


Vary 


H,S 


Vary Vary 


Vary 


Dry deposition velocity^ 






H2 







CH^ 









1.0 X 10~^ (Assumed) 




CO 


1.0 X 10~* em s~^ (Kliarecha et al. 2005) 




CHjO 


0.1 em s~l (Wagner et al. 2002) 




CO2 


1.0 X 10"* em s^^ (Archer 2010) 




O2 







O3 


0.4 em s~^ (Hauglustaine et al. 1994) 




H2O2 


0.5 cm s""*^ (Hauglustaine et al. 1994) 




HjS 


0.015 cm s-1 (Sehmel 1980) 




SO2 


1.0 cm (Sehmel 1980) 






0.2 cm (Sehmel 1980) 




H,SO^(A) 


0.2 cm (Sehmel 1980) 





"Rainout rates for Hg, CO, CH4, C2Hg, and O2 arc generally assumed to be zero to simulate an ocean surface 
saturated witli tliesc gases on an abiotic exoplanet. 

''The volcanic gas emission rates from the planetary surface are assigned for each model scenario. H2O 
emission is not explicitly considered because the surface has a large water reservoir, i.e., an ocean. 

'^We here list the dry deposition velocities (with references) for emitted gases and their major photochemical 
byproducts, and dry deposition velocities that are important for the mass and redox balance of the atmo- 
sphere. Dry deposition velocities are assumed to be identical for the three scenarios. C2Hg dry deposition 
velocity is assumed to take into account the loss of carbon due to organic haze formation and deposition. 
The CO2 dry deposition velocity is assumed such that the steady-state mixing ratio of CO2 in H2 and N2 
atmospheres is in the order of 100 ppm. 
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required all models to balance mass flux within 10^^ for convergence, and typically our 
models balance mass flux to 10~^. We have also explicitly checked the redox (i.e., hydrogen 
budget) balance (see our deflnition of the redox number, flux, and balance in Paper I), and 
required the models to balance redox flux to 10~^. 

One of the most signiflcant controlling factors that determine the steady-state compo- 
sition of atmospheres on terrestrial exoplanets is the dry deposition velocities of emitted 
gases and their major photochemical byproducts. Of particular importance in this paper is 
the dry deposition velocities of and SO2, which could vary by orders of magnitude (see 
Table [1] for the flducial values of key dry deposition velocities). The dry deposition velocities 
depend on the properties of the lower atmosphere and the surface. For example, in a model 
of the early cold Martian atmosphere, the deposition velocities are assumed to be reduced 
by an artificial factors of up to 1000 compared with those in warm current Earth (e.g., Tian 
et al. 2010), in order to account for less efficient deposition at a lower temperature. For an- 
other example, in modeling the early Martian atmosphere, because the putative Mars ocean 
is believed to be saturated with dissolved SO2 and other sulfur species, SO2 deposition is 
assumed to be balanced by the an equivalent return flux from the ocean (e.g. Halevy et 
al. 2007; Tian et al. 2010) and then Vy)ep of SO2 is assumed to be zero. In this paper, 
we explore the effect of varying the dry deposition velocities of H2S and SO2, which can be 
scaled down by a factor of 100 if the surface is saturated with sulflde or sulflte. 

Another important factor that strongly influences the atmospheric sulfur chemistry is 
the formation and sedimentation of aerosols. Photochemically produced H2SO4 and Sg may 
condense to form aerosols if their concentrations exceed their vapor saturation concentra- 
tions. The parameterization of condensation and sedimentation of aerosols in our photo- 
chemistry model is described in detail in Paper I. Saturation vapor pressure of H2SO4 is 
taken as recommended by Seinfeld & Pandis (2006) for atmospheric modeling, with a valid- 
ity temperature range of 150 - 360 K. Sg is the stable form of elemental sulfur because the 
Sg molecule has a crown-shape ring structure that puts least strain on the S-S bond among 
sulfur allotropes and the crown structure allows for considerable cross-ring interaction be- 
tween nonbonded atoms (Meyer 1976). The saturation pressure of Sg is then taken as the 
total sulfur saturation pressure against liquid sulfur at T > 392 K and solid (monoclinic) 
sulfur at T < 392 K tabulated by Lyons (2008). We consider the average aerosol particle 
diameter, a key parameter that determines the aerosols' dynamical and optical properties, to 
be a free parameter. On Earth, the ambient aerosol size distribution is dominated by several 
modes corresponding to different sources. The "condensation submode" , formed from vapor 
condensation and coagulation, has an average diameter of ~ 0.4 fim (Seinfeld & Pandis, 
2006). On Venus, the Mode 1 particles with an an average diameter of ~ 1.0 fim dominate 
the upper cloud (e.g., Carlson et al. 1993). On Titan, the photochemical aerosols in the 
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stratosphere have mean diameters in the range of 0.1 - 1 //m (Rages et al. 1983). We treat 
the particle diameter a free parameter, and explore the effects of varying the particle diam- 
eter from 0.1 to 10 /xm. The dry deposition velocity of aerosols is assumed to be 0.2 cm s~^, 
a sensible deposition velocity of particles having diameters between 0.1 and 1 /xm (Sehmel 
1980; Seinfeld & Pandis 2006). 

Before leaving this section we provide the physical rationale of specifying temperature- 
pressure profiles and eddy diffusion coefficients for our atmospheric photochemistry models. 

We modeled the atmospheric composition from the 1-bar pressure level up to the alti- 
tudes of about 10 scale heights. We chose appropriate vertical resolution for each scenario so 
that there are 4 layers per scale height. We adopted a temperature profile for our atmosphere 
models, without considering feedback on temperature of the atmospheric composition. The 
surface temperature is assumed to be 288 K. The semi-major axis of a terrestrial exoplanet 
around a Sun-like star implied by this surface temperature is 1.6 AU, 1.0 AU, and 1.3 AU, for 
H2, N2, and CO2 dominated atmosphere, estimated based on a similar procedure as Kasting 
et al. (1993). The temperature profiles are then assume to follow appropriate dry adiabatic 
lapse rate (i.e., the convective layer) until 160 K (Hg atmosphere), 200 K (N2 atmosphere), 
and 175K (CO2 atmosphere) and to be constant above (i.e., the radiative layer). The adopted 
temperature profiles arc consistent with significant greenhouse effects in the convective layer 
and no additional heating above the convective layer for habitable exoplanets. We did not 
consider the climate feedback of SO2 or sulfur aerosols in the atmosphere, which could be 
important to determine the surface temperature (e.g. Halcvy ct al. 2007; Tian et al. 2010). 
While not ideal, these temperature profiles yield the same results discussed in the below as 
temperature profiles varied by several tens of K. The precise temperature- pressure struc- 
ture of the atmosphere is less important than photochemistry for the investigation of sulfur 
chemistry because the most important photolysis and chemical reactions are not signifi- 
cantly affected by minor deviations in the temperature profile. We found that variation of 
temperature profiles by a few tens of K has minor impact on the atmospheric composition. 

Vertical transport of gases in the atmosphere is parameterized by eddy diffusion, and 
the coefficients are assumed to be those of Earth's atmosphere scaled by the atmospheric 
scale height to account for H2, N2, and CO2 being the dominant species. The current Earth's 
eddy diffusion coefficient profile has been derived from mixing ratio profiles of several long- 
lived gases (Massie & Huntcn 1981; also shown in Figure 1 of Paper I). We use the empirical 
eddy diffusion coefficient profile for current Earth as a template, and scale the coefficient 
inversely with the mean molecular mass for H2, N2, and CO2 dominated atmospheres. The 
justification for such scaling is that the eddy diffusion coefficient is proportional to the 
mixing length which is in turn a portion of the atmospheric scale height (e.g.. Smith 1998). 
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The pressure surface to pressure surface projection also ensures that the eddy diffusion 
coefficient profile features a eddy diffusion minimum near the tropopause for atmospheres 
with different mean molecular masses. Our approach to parameterize vertical transport is 
of course an approximation. An accurate representation of vertical transport would likely 
involve circulation on the global scale (e.g., Holton 1986). We explore the effect of eddy 
diffusion coefficients ranging one or two orders of magnitude from the nominal value in a 
sensitivity study (see section 4.1). 



2.2. Radiative Transfer Model 

We compute synthetic spectra of the modeled exoplanet's atmospheric transmission, 
reflection and thermal emission with a line-by-line radiative transfer code (Scagcr & Sasselov, 
2000; Seager et al. 2000; Miller-Ricci et al. 2009; Madhusudhan & Seager 2009). Opacities 
are based on molecular absorption with cross sections computed based from the HITRAN 
2008 database (Rothman et al. 2009), molecular coUision-induced absorption when necessary 
(e.g., Borysow 2002), Rayleigh scattering, and aerosol extinction are computed based on the 
Mie theory (e.g.. Van de Hulst 1981). The transmission is computed for each wavelength 
by integrating the optical depth along the limb path, as outlined in Seager & Sasselov 
(2000). The reflected stellar light and the planetary thermal emission are computed by the 
5-Eddington 2-stream method (Toon et al. 1989). We used the refractive index of Sg aerosols 
from Tian et al. (2010) for the UV and visible wavelengths and from Sasson et al. (1985) 
for infrared (IR) wavelengths. We used the refractive index of H2SO4 aerosols (assumed to 
be the same as 75% sulfuric acid solution) from Palmer & William (1975) for UV to IR 
wavelengths, and Jones (1976) for far IR wavelengths. 

The particle size distribution of aerosols controls their optical properties. We adopt the 
lognormal distribution as 



dN Nt 

-TTT = , exp 

dD y/2^D In a 



(In D - In Do) 
21nV 



(1) 



where dN is the number of particles per volume in the diameter bin dD, Nf is the total 
number density of particles, Dq is the median diameter of the particles, and a is the particle 
size dispersion (dcflned as the ratio of the diameter below which 84.1% of the particles lie 
to the median diameter). The lognormal distribution is a reasonable assumption because it 
provides a good fit to the particle size distribution measured in Earth's atmosphere (e.g., 
Seinfeld & Pandis 2006) , and a sensible particle size dispersion parameter for photochemically 
produced aerosols is in the range of 1.5 ~ 2.0 (Seinfeld & Pandis 2006). What is important 
for the radiative transfer model and the photochemistry model is the surface area mean 
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diameter Ds and the volume mean diameter Dy, respectively. The mean diameters are 
related to the median diameter as 



We use the surface area mean diameter Ds (referred to as "mean diameter" in the following) 
as the free parameter for specifying a particle size distribution, as it is relevant to the 
radiative properties of the particle population. The volume mean diameter Dy is useful in the 
conversion from mass concentration of the condensed phase (computed in the photochemistry 
model) to the number of aerosol particles for radiative transfer computation. The extinction 
cross sections of H2SO4 and Sg molecules in aerosols for various mean particular diameters 
are shown in Figure [T] 

Elemental sulfur aerosols and sulfuric acid aerosols have different optical properties at 
the visible and infrared (IR) wavelengths. In the visible, Sg aerosols have a larger cross 
section than H2SO4 aerosols (Figure [1]). For wavelengths less than 400 nm Sg aerosols 
are both reflective and absorptive. In the infrared, the cross section of Sg aerosols drops 
significantly with increasing wavelength unless the mean diameter is in the order of 10 /im. 
In contrast, H2SO4 aerosols have an enhancement of absorption at the MIR wavelengths 
(5-10 /xm) for all particle sizes (Figured]). 



We now briefly describe the most important processes of sulfur chemistry that occur 
in reducing and oxidizing atmospheres on rocky exoplanets. The primary sulfur emission 
from the planetary surface would be SO2 and H2S; they are either deposited back to the 
surface via dry or wet deposition, or converted into other forms of sulfur compounds in the 
atmosphere by photochemical reactions. One of the main purposes of this paper is to study 
the fate of sulfur gases emitted from the surface and their possible photochemical byproducts 
in the atmosphere. 

The fate of sulfur gases emitted from the surface is mainly controlled by the redox 
power of the atmosphere - the ability to reduce or oxidize a gas in the atmosphere. The 
redox power, in turn, is controlled by both the main component in the atmosphere (e.g., 
H2, N2, and CO2) and the surface emission and deposition of trace gases (i.e., H2, CH4, and 
H2S), as shown in Table [21 In the extreme cases of the atmospheric redox state, i.e., the H2- 
dominated atmospheres and the 02-rich atmospheres, the atmospheric redox power is surely 




(2) 




(3) 
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reducing or oxidizing, regardless of the nature of surface emission or deposition. However, 
for an intermediate redox state, the atmosphere would be composed of redox-neutral species 
like N2 and COg, and the redox power of the atmosphere can be mainly controlled by the 
emission and the deposition fluxes of trace gases from the surface. The higher the emission 
of reducing gases is, the more reducing the atmosphere becomes. 

We already know that sulfur gas emissions are effectively oxidized into sulfate, the 
most oxidized form of sulfur, in the oxic atmospheres such as the Earth's (e.g., Seinfeld 
& Pandis 2006). In anoxic atmospheres, which include the reduced atmospheres and the 
oxidized atmospheres, previous studies have shown that both elemental sulfur and sulfate 
could be formed (Kasting 1990; Pavlov & Kasting 2002; Zahnle et al. 2006; Hu et al. 2012). 
For this paper, we use dominated atmospheres as the representative cases for reducing 
atmospheres, and we use Ng, and CO2 dominated atmospheres as the representative cases 
for oxidized atmospheres that could be both reducing and oxidizing. We now describe the 
key sulfur chemistry processes in these atmospheres as the follows. 

Table 2: Redox power of atmospheres on rocky exoplanets. 



Type 


Main Component 


Redox Power 


Main Reactive Species 


Solar-System analogs 


Note 


Reduced 


H^, CO 


Reducing 


H 


None 




Oxidized 


N2, CO2 


Weakly reducing 
Weakly oxidizing 


H, OH, 
H, OH, 


None 

Mars, Venus 


The redox power is mainly controlled by the 
surface emission of trace gases (H^, CH^, H^S). 


Oxic 


0. 


Highly oxidizing 


OH, 


Earth 





3.1. Reducing Ha-Dominated Atmospheres 

Both H2S and SO2 emitted from the surface are efficiently converted into elemental 
sulfur in reducing H2 atmospheres. The major chemical pathways for sulfur compounds 
in the H2-dominated atmosphere and the results of photochemistry model simulations are 
shown in Figure |21 Atomic hydrogen produced from photodissociation of water vapor and 
H2S itself is the key reactive species that converts H2S and SO2 into elemental sulfur. The 
primary chemical loss for H2S in the atmosphere is via 

H2S + hz/ — ^ HS + H, (CI) 

and 

H2S + H HS + H2. (C2) 

The HS produced can then react with H again or with itself to produce elemental sulfur. HS 
can also react with S to produce S2. The primary chemical loss for SO2 in the atmosphere is 
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photodissociation that produces SO. SO can be either photodissociated to elemental sulfur, 
or be further reduced to HS via HSO by H or CHO and then converted to elemental sulfur 
(see Figure [2]) ■ 

The S and S2 molecules produced in the atmosphere will polymerize to form Sg, and Sg 
will condense to form aerosols if it is saturated in the atmosphere. Due to its ring structure, 
Sg is stable against photodissociation. Sg is a strong UV absorber (Kasting 1990). Therefore, 
Sg aerosols, if produced in the atmosphere, can effectively shield UV photons so that H2S 
and SO2 may accumulate beneath the aerosol layer (see the case for a sulfur emission rate 
300 times higher than Earth's current volcanic sulfur emission rate shown in Figure |2]). 

The primary source of atomic hydrogen is the photolysis of H2O, which occurs above 
the altitudes of ~ 10'^ Pa pressure level. The atomic hydrogen can be then transported by 
eddy diffusion to the pressure level of ~ 10^ Pa to facilitate the removal of HgS and SO2 and 
the production of elemental sulfur. Additional numerical simulations show that an increase 
of the eddy diffusion coefficient by one order of magnitude can increase the yield of elemental 
sulfur by about 20%, because the transport of atomic hydrogen becomes more efficient. A 
secondary source of atomic hydrogen is photolysis of H2S (reaction IC ip . This secondary 
source for atomic hydrogen is particularly important when the host star is a quiet M dwarf, 
because a quiet M dwarf produces few photons that could dissociate watecl. For planets 
around quiet M dwarfs the photolysis of H2S could be the main source of atomic hydrogen 
in their atmospheres. Additional numerical simulations show that in the habitable zone of a 
quiet M dwarf having an effective temperature of 3100 K, H2S photolysis alone can produce 
enough atomic hydrogen to drive the formation of elemental sulfur in the atmosphere. 

We here comment on the uncertainty of photochemistry models regarding the yield of 
Sg. In our model, we have assumed polymerization of elemental sulfur proceeds via 



S + S— ^Ss, (C3) 

S + S2-^S3, (C4) 

S + S3-^S4, (C5) 

S2 + S2-^S4, (C6) 

S4 + S4 V Sg . (C 7) 



We also include photodissociation for S2, S3, and S4. However, the reaction rates of sulfur 
polymerization (reactions IC 31 - IC 71) have not been well established by laboratory studies. 



■^Water is principally dissociated by photons in the 150 - 200 nm wavelength range, whereas H2S is 
principally dissociated by photons in the 200 - 260 nm wavelength range. 
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and previous authors have adopted different rate constants for these reactions. In particular, 
Kasting (1990) and Pavlov & Kasting (2002) have used 3-order-of-magnitude lower rates for 
reactions f IC 3I - [(T4|) and 1-order-of-magnitude lower rates for reactions f IC 5I - [UT|) . compared 
with Moses et al. (2002). In this work, we have adopted the reaction rates of Moses et al. 
(2002) for elemental sulfur reactions. Our sensitivity tests show that adopting the reaction 
rates of Kasting (1990) would result in about 3 to 10 times less Sg. We have chosen a 
higher sulfur polymerization rates for nominal models because: (1) the chemical pathways of 
reactions (IC 3I - [(T7|) are probably not complete and there may be other pathways to form Sg; 
(2) Sg, Sg, and S4 may condense as suggested by Lyons (2008) and the polymerization may 
still proceed in the condensed phase to Sg. Experimental studies are encouraged to settle 
this important uncertainty. 

3.2. Oxidized N2 and COa-Dominated Atmospheres 

HgS and SOg gases emitted from the surface can be converted into both elemental 
sulfur (Sg) and sulfuric acid (II2SO4) in the oxidized (but anoxic) atmospheres such as Ng 
and CO2 dominated atmospheres. The major chemical pathways that lead to formation of 
both elemental sulfur and sulfuric acid, and the results of photochemistry model simulations 
are shown in Figure El The production of elemental sulfur aerosols involves UV photons and 
atomic hydrogen, as does in reducing Hg atmospheres; whereas the production of sulfuric 
acid requires oxidizing species, notably OH and Og. These reactive species, either reducing 
or oxidizing, are produced from photolysis of water and COg. In particular, the source of OH, 
responsible for converting SOg to sulfuric acid in the atmosphere, is the photodissociation of 
HgO. As a result, the amount of UV photons that are capable of dissociating water controls 
the yield of H2SO4. For example in the habitable zone of a quiet M dwarf the yield of H2SO4 
is much reduced compared with solar-like stars by at least one order of magnitude. 

The photochemically produced Sg and H2SO4 may condense to form aerosols in the 
atmosphere if saturated. As a result, aerosols in the atmospheres provide a UV shield that 
enables the accumulation of HgS and SO2 beneath the layer of aerosols. In particular for 
an Earth-sized planet in the habitable zone of a Sun-like star, when the surface emission 
rate is more than two-orders-of-magnitude higher than the current Earth's volcanic sulfur 
emission rate, photochemical aerosols in the atmosphere lead to substantial UV shielding for 
accumulation of HgS and SO2, as shown in Figure El We find that only when sulfur emission 
is highly elevated with respect to current Earth could HgS or SO2 accumulate to the order 
of parts per million mixing ratio in the Ng and CO2 atmospheres. 

The relative yield between elemental sulfur and sulfuric acid is controlled by the redox 
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power of the atmosphere. In general, more sulfuric acid aerosols and less elemental sulfur 
aerosols are anticipated in a more oxidizing atmosphere. The redox power of the atmosphere, 
in turn, is determined by both the main constituent and the reducing gas emission. COg 
dominated atmospheres are mosre oxidizing than dominated atmospheres as photodis- 
sociation of COg leads to atomic oxygen. Therefore the primary sulfur emission is more 
likely to be converted to sulfuric acid in CO2 dominated atmospheres than in dominated 
atmospheres (see Figure |3]). Surface emission of reducing gases, including Hg, CH4, CO, and 
H2S, alters the redox budget of the atmosphere and therefore increases the relative yield of 
elemental sulfur versus sulfuric acid aerosols (e.g., Zahnle et al 2006). As shown in Figure 
|3l when the sulfur emission rate increases, both and CO2 atmospheres become more and 
more reducing (because H2S is reducing), which results in a dramatic increase of elemental 
sulfur production in the atmosphere. Furthermore, the HgS/SOg ratio in the sulfur emission 
affects its contribution to the redox power of the atmosphere and then the relative abun- 
dances of the two types of aerosols in the atmosphere significantly. As a result of the increase 
in the H2S/SO2 emission ratio, the amount of Sg aerosol in the atmosphere increases, and the 
amount of H2SO4 in the atmosphere decreases dramatically (see Figure H]). For an Earth-like 
planet having an N2 atmosphere, if the H2S/SO2 emission ratio is less than 0.1 (as is the case 
for current Earth; Holland 2002), the dominant type of aerosols in the atmosphere is sulfate; 
whereas elemental sulfur aerosols become the dominant type if the H2S/SO2 emission ratio 
is larger than 1. 



4. Results 

4.1. Optically Thick Aerosols from Sulfur Emission 

The main finding is that on terrestrial exoplanets having atmospheres ranging from 
reducing to oxidizing, the primary sulfur emission from the surface (e.g., H2S and SO2) is 
chemically short-lived. The sulfur emission leads to photochemical formation of elemental 
sulfur (Sg) and sulfuric acid (H2SO4), which would condense to form aerosols if saturated in 
the atmosphere. In reducing atmospheres (e.g., H2 atmospheres), Sg aerosols are photochem- 
ically formed based on H2S and SO2 emission; and in oxidized atmospheres (e.g., N2 and 
CO2 atmospheres), both Sg and H2SO4 aerosols may be formed (see Figure E]). In general, 
the higher the surface sulfur emission, the more aerosols exist in the atmosphere (see Figure 

ED. 

As a result of photochemical production of elemental sulfur and sulfuric acid, terrestrial 
exoplanets with a habitable surface temperature (e.g., 270 ~ 320 K) and substantial sulfur 
emission from the surface are likely to have hazy atmospheres. In this paper, we use "hazy" 
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to describe an atmosphere that has significant aerosol opacities at visible wavelengths (e.g., 
500 nm). We find that even with an Earth- like surface sulfur emission, 1-bar H2 dominated 
atmospheres on habitable rocky exoplanets are hazy with Sg aerosols (see Figure E]). We also 
find that if the sulfur emission rate is 30 ~ 300 times more than the Earth's current volcanic 
sulfur emission rate, photochemical Sg aerosols become optically thick at visible wavelengths 
in oxidized atmospheres including and CO2 atmospheres (see Figure [5]). 

The key parameters that determine the aerosol opacity in the atmosphere are the surface 
sulfur emission rate, the dry deposition velocity, and the aerosol particle size. First, a higher 
surface sulfur emission rate leads to more sulfur and sulfate aerosols in anoxic atmospheres 
(e.g.. Figure [5] and Figure [6]). Second, larger dry deposition velocities of HjS and SO2 
cause more rapid removal of these sulfur compounds from the atmosphere, which reduces 
the chance of converting them into condensable molecules (i.e., Sg and H2SO4). Therefore, 
larger dry deposition velocities of H2S and SO2 result in lower aerosol loading and aerosol 
opacities in the atmospheres, as shown in Figure [71 Third, we find that the particle size has 
only secondary effects on the chemical composition of the atmosphere (i.e., by increasing the 
penetration of ultraviolet radiation), but has a primary effect on the aerosol optical depth. 
For mean particle diameter varying in the range of 0.1 ~ 1 fim (i.e., typical particle sizes 
of photochemical aerosols on Earth (e.g., Seinfeld & Pandis 2006) and Titan (e.g.. Rages 
et al. 1983)), we do not see a notable variation in the yield of elemental sulfur, but we 
see an enhancement of H2SO4 production with large particles (Figure [6]). Even with the 
same aerosol abundances, however, micron-sized particles cause lower opacities at the visible 
wavelengths and higher opacities in MIR compared with submicron-sized particles (Figure 



We capture the effects of the three key parameters on the aerosol opacity in anoxic 
atmospheres on terrestrial exoplanets by fitting the following power-law formula, i.e.. 



where r is the vertical optical depth due to aerosols at 1 bar, $(S) is the total sulfur emission 
rate, Vdep/Hdep (Earth) is the dry deposition velocities of H2S and SO2 with respect to 
current Earth values, dp is the mean particle diameter of aerosols, a, b, and c are positive 
numbers, and C is a constant that covers other uncertainties. We have fit the empirical 
relation (jl]) through an extensive parameter exploration using photochemistry models (see 
Figure [5] - [7| for examples) and determined the values of C, a, b and c for H2 dominated 
reducing atmospheres and for N2 and CO2 dominated oxidized atmospheres. We summarize 
graphically the parameter regime in which sulfur emission leads to a hazy atmosphere in 
Figure [3 Here we use rsoonm and ry.s^m as the representatives for aerosol opacities at visible 



ED. 
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wavelengths and MIR wavelengths; due to the complex nature of the extinction cross sections 
of aerosol particles (Figured]), it is not practical to fold the full wavelength dependency into 
the empirical formula. Also, we find that it is always true that rsoonm ~ for mean particle 
diameter in the order of 10 /im and tj^^^^ ~ for mean particle diameter in the order of 0.1 
/im. 

For H2 atmospheres, and mean particle diameter d-p in the the range of 0.1 ~ 1 /xm, 

X -1 on^ ^(S) VY ^DEP V°V y ... 



and for dp in the the range of 1 ~ 10 fim 



-m ^^^^ ^°^v ^^^p ^"°^v (f^^ 

Vl0iicm-2s-iy' V^DEp(Earth)y' [l.O fimj ' 

For N2 and CO2 atmospheres, and dp in the the range of 0.1 and 1 /im, 

($(S) \ 0.7 / ^ \ — 0.3 / ^ \ — 0.7 

1011 cm-2 s-i J (vDEP(Earth)J (oTT^ j ' 



and for dp in the the range of 1 ~ 10 /im 

r7.5Mm = 0.01 ~ 0.1 



^/o\ \ 0.8 / \ -0.5 / T \ -1.6 

$(8) \ / V^DEP \ / «P 



10" cm-2 s-1 y V^DEP (Earth) / VO-1 yum 

The constant C in equations ([5] -[8]) spans about one order of magnitude, which covers 
the variation of the following model inputs: 

• The H2S/SO2 ratio in the surface sulfur emission, ranging from 0.01 to 10; 

• Temperature profiles deviating from the adopted temperature profile by ±30 K that 
controls the mixing ratio of water vapor in the atmosphere by the cold trap; 

• Stellar ultraviolet radiation received by the planet, ranging from the habitable zone of 
solar-like stars to the habitable zone of quiet M dwarfs with an effective temperature 
of 3100 K; 

• Eddy diffusion coefficients ranging from 0.1 to 100 times the values of Earth's atmo- 
sphere; 



Sulfur polymerization reaction rates (reactions IC 31 - [CT|) ranging by one order of mag- 
nitude. 
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To summarize, we find that the emission of H2S and SO2 from the surface is readily 
converted into sulfur (Sg) and sulfate (H2SO4) in anoxic atmospheres of terrestrial exoplanets. 
The photochemical sulfur and sulfate would condense to form aerosols if saturated in the 
atmosphere, which is likely to occur on a planet in the habitable zone of either a Sun-like star 
or a quiet M star. The aerosol layer is optically thick at the visible and NIR wavelengths 
if the surface sulfur emission is comparable to Earth's volcanic sulfur emission in the 
atmosphere, and more than 30 ~ 300 times of the Earth's volcanic sulfur emission in other 
anoxic atmospheres, depending on the dry deposition velocities of sulfur compounds and 
particle size of the aerosols. 



4.2. Spectral Features of SO2, H2S, and Sg and H2SO4 Aerosols 

The sulfur emission from surface shapes the spectra of terrestrial exoplanets at the visible 
and NIR wavelengths, mostly through the photochemical formation of Sg and H2SO4 aerosols. 
We use the model outputs from the photochemistry models to compute the transmission, 
refiection, and thermal emission spectra of a terrestrial exoplanet with various levels of 
sulfur emission, and show examples of the computed spectra in Figure |9l We see that 
submicron-sized Sg aerosols dominate the transmission and refiection spectra at wavelengths 
from visible up to 3 /im, if the sulfur emission is more than about two orders of magnitude 
higher than Earth's volcanic sulfur emission. In general, an atmosphere with high sulfur 
emission and therefore high aerosol loading generally exhibits a fiat transmission spectrum 
(the H2O features at NIR muted), and a high visible albedo (see Figure [9]). Notably, Sg 
aerosols are purely refiective at 500 nm but absorptive at 300 nm. The absorption edge of Sg 
aerosols in 300 - 400 nm is evident in the refiection spectra for planets with enhanced sulfur 
emission (Figure [H]), which is a potential diagnostic feature for Sg aerosols. 

Although opaque at visible wavelengths, the atmospheres with enhanced sulfur emission 
are likely to be transparent in the MIR wavelengths (A > 5 fim). The spectral features of 
aerosols depend on their particle sizes, so we now consider two possibilities: if the particles 
are submicron-sized, the aerosol molecules have negligible cross sections at MIR (see Figure 
[1]); or if the particles are micron-sized, the falling velocity of aerosol particles is large enough 
to rapidly remove aerosols from the atmosphere, as implied by equation ([6]) and equation ([8]) 
that are applicable for micron-sized particles. Therefore in both cases the aerosol opacities 
at MIR are minimal even for very high sulfur emission rates (see Figure M for examples 
of N2 atmospheres, and atmospheres are qualitatively similar). The only exception, in 
which aerosols indeed affect MIR spectra, is the case of abundant H2SO4 aerosols. The main 
spectral effect of H2SO4 aerosols is absorption at MIR wavelengths (5 ~ 10 /xm; Figure [9]). 
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However, the column- average mixing ratio of H2SO4 needs to be larger than 0.1 ppm in order 
to produce significant aerosol absorption at MIR. We find with numerical exploration that 
such a high abundance of H2SO4 aerosols is only possible in highly oxidizing CO2 dominated 
atmospheres without reducing gas emission (see Paper I for an example of such atmospheres). 
With reducing gas emission (i.e., and CH4), it is unhkely that H2SO4 mixing ratio exceeds 
0.01 ppm in anoxic atmospheres for a wide range of sulfur emission rates (see Figure |5]). In 
summary we expect the spectral effects of Sg and H2SO4 aerosols to be minimal at MIR for 
most cases. 

We now turn to consider the direct spectral features of H2S and SO2. It has been 
previously proposed that H2S and SO2 can be detectable on terrestrial exoplanets by their 
spectral features (Kaltenegger & Sasselov 2010). However, our photochemistry models show 
that both HgS and SO2 are chemically short-lived in the atmospheres, which implies that 
that substantial surface emission is required to maintain a detectable level of either H2S or 
SO2 in the atmosphere. SO2 has diagnostic absorption features at 7.5 fim and 20 yum (see 
Figure [9]). For these features to be detectable the mixing ratio of SO2 needs to be larger 
than 0.1 ppm, which corresponds to sulfur emission rates 1000 times more than current 
Earth's sulfur emission rates for H2, N2, and CO2 atmospheres (see Figure E]). The spectral 
feature of HgS is the pseudo-continuum absorption at wavelengths longer than 30 fim, which 
coincides with the rotational bands of H2O. We find that the only scenario in which HgS may 
be directly detected is the case with extremely high sulfur emission rates (i.e., 3000 times 
higher than the current Earth's sulfur emission rate) on a highly desiccated planet without 
liquid water ocean so that there is no water vapor contamination. We therefore conclude that 
direct detection of H2S and SO2 is tricky: they are chemically short-lived so that extremely 
large surface emission is required for a detectable mixing ratio in the atmosphere, and their 
spectral features may be contaminated by other gases in the atmosphere. 

Finally, we suggest that the emission of sulfur compounds might be indirectly inferred 
by detecting sulfur and sulfate aerosols. Our numerical exploration reveals a monotonic 
relationship between the abundance of aerosols in the atmosphere and the emission rates 
of sulfur compounds (see Figure [5]), and the composition of aerosols is correlated with the 
H2S/SO2 ratio of the surface emission (see Figure Hj). A combination of featureless low 
atmospheric transmission (large planet radius viewed in transits) and high planetary albedo 
(large planetary fiux at the visible wavelengths viewed in occultations) may establish the 
existence of aerosols in the atmosphere. In particular, elemental sulfur (Sg) aerosols are 
absorptive at wavelengths shorter than 400 nm and therefore might be identified by the 
absorption edge (see Figure [2]). Sulfate aerosols (H2SO4), if abundant in the atmosphere, 
lead to absorption features at the MIR wavelengths (A ~ 5 - 10 fim). However, none of 
these features are uniquely diagnostic of certain types of aerosols. The identification of 
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aerosol composition, therefore, is by no means straightforward. We learn from the Solar 
System exploration that the discriminating piece of information for aerosol identification 
comes from polarization of reflected stellar light. Historically, the bright clouds on Venus 
were identified to be mainly composed of H2SO4 droplets after the phase curve of the planet 
in polarized light had been observed (e.g.. Young 1973; Hansen & Hovenier 1974). We 
therefore postulate that aerosol identification on terrestrial exoplanets and the inference of 
surface sulfur emission might require observation of polarized reflected light as a function of 
planetary illumination phase. 



Biosignatures are gases in an exoplanet's atmosphere produced by life. In order to 
confirm a certain gas to be a plausible biosignature, one has to verify that the gas can 
accumulate in the atmosphere and that the amount of gas detected cannot be produced 
through plausible abiotic processes. It has been proposed and widely discussed that O2 
(and its photolytic product O3), N2O, CH4 and CH3CI are exoplanet biosignatures (e.g., 
Sagan et al. 1993; Des Marais et al. 2002; Segura et al. 2005; Segura et al. 2007; Hu 
et al. 2012). Organosulfur compounds, such as dimethyl sulfide (DMS) and methanethiol 
(CH3SH), have also been suggested to be biosignatures of the early Earth (Pilcher 2003) and 
anoxic exoplanets (Domagal-Goldman et al. 2011). 

HgS can be produced from several metabolic origins on Earth, and so is a candidate 
biosignature gas. Life on Earth can produce H2S through sulfate reduction (when the en- 
vironment is reduced) and sulfur disproportionation. Microorganisms can disproportionate 
sulfur compounds of intermediate oxidation states, including thiosulfate, sulfite, and elemen- 
tal sulfur, into H2S and sulfate (Finster 2008). For example, disproportionation of sulfite in 
the ocean is described by 



in which the Gibbs free energy released is 58.9 kJ mol~^ sulfite. The sulfite reducers, in- 
cluding microorganisms in genus Desulfovibrio and Desulfocapsa, extract energy from the 
disproportionation (Kramer and Cypionka, 1989). 

The effect of biotic HgS production is the increase of the HgS/SOg ratio of the surface 
sulfur emission. If the H2S/SO2 ratio in the volcanic sulfur emission is low (i.e., less than 
0.1), sulfur disproportionation and sulfate reduction by life may increase the H2S/SO2 ratio 
significantly, which may lead to a change in the redox input to the atmosphere and therefore 
the dominant aerosol species in the atmosphere, as suggested by Figure HI Specifically, for a 



5. Discussion: Can H2S be a Biosignature Gas? 




+ HS 



(C8) 
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habitable terrestrial exoplanet having a weakly oxidizing atmosphere, biotic production 
of in excess of the geological HgS emission could result in a higher amount of Sg aerosols 
and a much lower amount of H2SO4 aerosols in the atmosphere compared with a planet 
without life. Although it is currently not possible to distinguish different types of aerosols, 
HgS could be a biosignature gas in the long term. 

The geological production of H2S, and the consequent risk of a false positive mis- 
identification of geological HgS for biological H2S, will be a major obstacle of confirming 
H2S to be a biosignature gas. Sulfur is believed to be present in the mantle of all terrestrial 
planets, and what determines the H2S/SO2 ratio in the volcanic outgassing is the oxygen 
fugacity of the upper mantle, temperature of the location where magma degassing happens, 
water content in the conduit of magma, and gas content dissolved in the magma (e.g., Hol- 
land 1984; Kasting et al. 1985; Holland 2002; Burgisser & Scaillet 2007). The current Earth 
volcanic emissions are oxidized, dominated by HgO, CO2 and SO2 with minor contributions 
of H2, CO and H2S. This volcanic gas composition is consistent with a magma buffered 
by the quartz-fayalite-magnetite (QFM) equilibrium, i.e., a relatively oxidized upper mantle 
(Holland 1984). As a global average, the volcanic H2S/SO2 emission ratio on Earth is 0.1 
(Holland 2002). However, if the mantle of a rocky exoplanet is much more reducing than 
that of the current Earth, significant geological source of H2S can be expected (Holland 1984; 
Kasting 1993). As a result, it would be very hard to rule out a geological contribution to 
the H2S emission flux by remote sensing. 

In summary, although H2S can be produced by energy-yielding metabolism, it is very 
unlikely to be a useful biosignature gas for three reasons. Firstly, H2S itself is unlikely to 
be detectable directly by remote sensing because of its weak spectral features, and their 
contamination by the spectral features of water. Secondly, HgS has a very short atmospheric 
lifetime, and so unrealistic emission rates are required to build up significant levels in any 
atmosphere. This second point could be overcome, in principle, by detecting Sg aerosols 
in an anoxic atmosphere, and discriminating them from H2SO4 aerosols. Discrimination 
between sulfur and sulfuric acid aerosols is not possible with current equipment, but may be 
possible in the future through analysis of refiected light. Thirdly, however, H2S suffers from 
a significant false positive risk, as geological sources can also produce H2S, and the ratio of 
H2S/SO2 in geological emissions depends on mantle chemistry, the physical structure of the 
outgassing events, and the extent of surface reprocessing of vented sulfur gases. To infer 
that life was generating H2S on a planet, this study shows that the observer would have to 
determine the Sg/H2S04 aerosol ratio and have knowledge of the geological outgassing ratio 
of H2S/SO2 and have knowledge of the surface chemistry that might modulate the primary 
outgassing rate. This seems an unreasonable requirement. 
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6. Conclusion 

We studied the effect of and SO2 surface emission on anoxic atmosplieres of ter- 
restrial exoplanets. Witli a newly established one-dimensional photochemistry model that 
treats all relevant chemical reactions and photochemical processes of O, H, C, and S bearing 
species, as well as formation and sedimentation of sulfur and sulfate aerosols, we find that 
H2S and SO2 gases emitted from surface are chemically short-lived in both reducing and 
oxidizing atmospheres. The sulfur emission results in photochemical production of elemen- 
tal sulfur (Sg) and sulfuric acid (H2SO4), which would condense to form aerosols if they 
are saturated in the atmosphere. For a planet in the habitable zone of a Sun-like star or 
a M star, Earth-like sulfur emission rates would result in optically thick aerosol layers in 
H2-dominated atmospheres; and a sulfur emission rate 2-orders-of-magnitude higher than 
the Earth's volcanic sulfur emission rate would result in optically thick aerosol layers in 
N2 and COs-dominated atmospheres. The composition of the photochemically produced 
aerosols mostly depends on the redox state of the atmosphere: Sg aerosols are formed in the 
reducing atmospheres (e.g., atmospheres), and both Sg and H2SO4 aerosols are formed 
in the oxidized atmospheres that could be both reducing and oxidizing (e.g., and CO2 
atmospheres). Based on extensive numerical simulations, we provide empirical formulae that 
show the dependency of the aerosol optical depth on the surface sulfur emission rates, the 
dry deposition velocities of sulfur compounds, and the aerosol particle sizes. 

We find that direct detection of H2S and SO2 is unlikely due to the rapid photochemical 
conversion from H2S and SO2 to elemental sulfur and sulfuric acid in atmospheres having a 
wide range of redox powers. For a terrestrial exoplanet with sulfur emitted from the surface 
at an enhanced rate, it is likely that at visible wavelengths the planet's atmosphere appears 
to be opaque due to the aerosol loading and that the planet has high visible albedo. However, 
for Earth-like planets with 1-bar atmospheres ranging from reducing to oxidizing, we find 
the effect of photochemical sulfur and/or sulfate aerosols in the MIR wavelengths is mini- 
mal, because micron-sized particles that interact with MIR photons have large gravitational 
setthng velocities and therefore short atmospheric hfetime. Finally, as the aerosol composi- 
tion is tightly related to the ratio of the H2S versus SO2 emission, although direct detection 
of H2S and SO2 by their spectral features is unlikely, their existence might be inferred by 
observing aerosol- related features in reflected light with future generation space telescopes. 
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Fig. 1. — Extinction cross sections (black lines) and scattering cross sections (orange lines) 
of H2SO4 and Sg per molecule in the condensed phase. The dotted, solid, and dashed lines 
are cross sections for the mean particle diameter of 0.1, 1, 10 /im, respectively. The size 
distribution of aerosol particles is assumed to be lognormal with a dispersion cr = 2. 
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Fig. 2. — Formation of elemental sulfur aerosols in reducing H2-dominated atmospheres 
on an Earth-sized rocky planet in the habitable zone of a Sun-like star. The upper panel 
schematically illustrates the chemical pathways from the primary sulfur emission (i.e., HgS 
and SO2) to elemental sulfur, in which sulfur compounds are located according to their 
oxidation states labeled on the top of the figure. The thin arrows show the major chemi- 
cal pathways in the atmosphere, and the thick arrows show the major surface-atmosphere 
interactions. The lower panel shows the results of photochemistry simulations, with total 
surface sulfur emission of 10^'^ (solid lines) and 10^^ (dashed lines) cm~^ s~^, i.e, 3 and 300 
times higher than the Earth's volcanic sulfur emission rate. The HgS/SOg ratio in the sulfur 
emission is 0.5 and the particle mean diameter is 0.1 fim. Other model parameters are tab- 
ulated in Table [1] UV photons and atomic hydrogen effectively convert H2S and SO2 into 
elemental sulfur, and elemental sulfur aerosols shield UV photons so that H2S and SO2 may 
accumulate below the aerosol layer if the sulfur emission is more than 300 times higher than 
the Earth's volcanic emission rate. 
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Fig. 3. — Formation of sulfuric acid aerosols in oxidized and COg dominated atmospheres 
on an Earth-sized rocky planet in the habitable zone of a Sun-like star. Similar to Figure 
[21 the upper panel schematically illustrates the chemical pathways from the primary sulfur 
emissions (i.e., H2S and SO2) to elemental sulfur and sulfuric acid. For double arrows the 
label above the arrow indicates the oxidizing agents, and the label below the arrow indicates 
the reducing agents. The lower two panel shows the results of photochemistry simulations for 
N2 and CO2 dominated atmospheres, respectively. The total surface sulfur emission is 10^*^ 
(solid lines) and 10^^ (dashed lines) cm~^ s~^, i.e, 3 and 300 times higher than the Earth's 
volcanic sulfur emission rate. The H2S/SO2 ratio in the sulfur emission is 0.5 and the particle 
mean diameter is 0.1 fim. Other model parameters are tabulated in Tabled] Both elemental 
sulfur aerosols and sulfuric acid aerosols are formed in the oxidized and anoxic atmospheres. 
The origins of the principle reducing agents (H and CHO) and the principle oxidizing agents 
(OH, O and O2) is photodissociation of H2O and CO2. The apparent depletion of SO2 at the 
pressure level of 10 - 100 Pa in the N2 atmosphere (the blue solid line in the middle panel) 
is due to the production of atomic hydrogen by methane photodissociation at this pressure 
level. 
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Fig. 4. — Correlation between the aerosol composition and the composition of sulfur emis- 
sions in the weakly oxidizing N2 atmosphere on an Earth-sized rocky planet in the habitable 
zone of a Sun- like star. The total surface emission rate is 10^^ cm~^ s~^, or 30 times the 
Earth's volcanic sulfur emission rate, and the particle mean diameter is 0.1 fim. Other model 
parameters are tabulated in Table [H As a larger fraction of surface sulfur emission is in the 
form of HgS, the amount of Sg aerosols in the atmosphere increases, and the amount of 
H2SO4 aerosols in the atmosphere decreases dramatically. 
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Fig. 5. — The relationship between the aerosol mixing ratios and aerosol opacities and 
the total sulfur emission rate. Column-integrated aerosol opacities at the 1-bar pressure 
level including both elemental sulfur aerosols and sulfuric acid aerosols at 500 nm (solid 
lines) and 7.5 fim (dashed lines) are shown in the lower panel. The planet is an Earth- 
sized rocky planet orbiting a Sun-like star, with reducing (Hg-dominated), weakly oxidizing 
(Ng-dominated), or highly oxidizing (COg-dominated) atmospheres. The aerosol particle 
mean diameter is assumed to be 0.1 /im, and the HgS/SOg ratio of the surface emission 
is 0.5. Other model parameters are tabulated in Table [H Sulfur emission 2-orders-of- 
magnitude larger than current Earth's volcanic sulfur emission (~ 3x 10^ S cm~^ s~^) leads to 
substantial aerosol opacities in the visible wavelengths in N2 and CO2 atmospheres, and sulfur 
emission comparable with current Earth's volcanic sulfur emission leads to substantial aerosol 
opacities in the visible wavelengths in atmospheres. The wiggle in the concentration of 
sulfuric acid aerosols reflects the competition between two effects: more sulfur is available to 
be converted into sulfuric acid as the sulfur emission increases, but the atmosphere becomes 
more reducing and less oxidizing as the sulfur emission increases. 
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Fig. 6. — Aerosol mixing ratios and optical depths at the surface (1-bar pressure level) at 
500 nm and 7.5 /im as a function of total sulfur emission rates, for an Earth-sized rocky 
planet orbiting a Sun-like star at 1 AU (black lines), a habitable planet around quiet M 
dwarf having effective temperature of 3100 K (red lines), and an Earth-sized rocky planet 
orbiting a Sun- like star at 1 AU with particle mean diameter of 1 fim (blue lines). The left 
column shows the case of weakly oxidizing N2 atmospheres, and the right column shows the 
case of reducing H2 atmospheres. The H2S/SO2 ratio in the sulfur emission is 0.5 and other 
model parameters are tabulated in Table [TJ We see that decreasing UV photon flux has little 
effect on the Sg formation, but results in a decrease of the amount of sulfuric acid aerosols, 
and therefore a decrease of MIR optical depth. We also see that particle diameter variation 
in 0.1 ~ 1 fim has little effect on the chemical composition, but for similar mass abundance 
a larger particle size results in a smaller optical depth in the visible wavelengths and a larger 
optical depth in the MIR wavelengths. 
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Fig. 7. — The relationship between the aerosol opacities (both Sg and H2SO4 aerosols) at 
500 nm and 7.5 fim in the H2 and atmospheres and the H2S and SO2 dry deposition 
velocities. Model parameters are shown in the figure and tabulated in Table [H 
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Fig. 8. — Formation of optically thick aerosols in atmospheres on rocky exoplanets in 
the habitable zone of their host star as a result of surface sulfur emission and deposition. 
The shaded areas are the parameter regime boundaries between a clear atmosphere and an 
optically thick atmosphere (defined as aerosol optical depth at 500-nm wavelength rsoonm > 
1), for reducing (Hg) and oxidized (Ng and COg) atmospheres. The upper-left corner of 
the parameter regime (small sulfur emission rates, large deposition velocities) leads to clear 
atmospheres; whereas the lower-right corner of parameter regime (large sulfur emission rates, 
small deposition velocities) leads to optically thick aerosols in the atmosphere composed of 
sulfur (Sg) and sulfate (H2SO4). The widths of the shaded boundary regime between clear 
atmospheres and optically thick atmospheres contain the uncertainties of: (1) the mean 
aerosol particle size ranging from 0.1 to 1 yum, (2) the HgS/SOg ratio of the sulfur emission 
ranging from 0.01 to 10 (i.e., more leads to thicker haze), (3) the spectral type of the host 
star ranging from G2V to M5, (4) the strength of vertical mixing in the atmosphere by eddy 
diffusion ranging from 0.1 to 10 times Earth's value, and (5) the surface temperature ranging 
from 270 to 320 K. Earth and Venus are shown for a reference in the Solar System: Earth's 
volcanic sulfur emission and HgS deposition velocity are plotted; and Venus' equivalent 
upward SO2 flux and SO2 deposition velocity at the altitude of 58 km are plotted (adapted 
from Krasnopolsky 2012). We note that the equivalent SO2 flux is a transfer rate across the 
58-km altitude, and does not imply a surface emission rate. The SO2 features at 7.5 fim 
and 20 fim requires a mixing ratio on the order of ppm to be spectrally significant, which 
corresponds to a sulfur emission flux of more than 10^^ cm^^ s~^ due to rapid photochemical 
removal of SOg in the atmosphere. 
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Fig. 9. — Transmission, reflection and thermal emission spectra of a terrestrial exoplanet 
with N2-dominated atmosphere with various surface sulfur emission up to 3000 times Earth's 
current volcanic emission (labeled in colors). The planet is an Earth-sized planet at the 1-AU 
orbit a Sun-like star, having surface temperature of 288 K. The H2S/SO2 emission ratio is 
0.5, the aerosol particle diameter is assumed to be 0.1 fim, and other model parameters are 
tabulated in Tabled! The cross sections of Sg and H2SO4 aerosols are shown in Figure [H 



